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Abstract. We have investigated the sliding of droplets made of solutions of Xanthan, a stiff rodlike polysaccharide 
exhibiting a non-newtonian behavior, notably characterized by a shear-rate dependence of the viscosity. These exper¬ 
imental results are quantitatively compared with those of newtonian fluids (water). The impact of the non-newtonian 
behavior on the sliding process was shown through the relation between the average dimensionless velocity (i.e. the 
Capillary number) and the dimensionless volume forces (i.e. the Bond number). To this aim, it is needed to define 
operative strategies to compute the Capillary number for the shear thinning fluids and compare with the corresponding 
newtonian case. Results from experiments were complemented with lattice Boltzmann numerical simulations of sliding 
droplets, aimed to disentangle the influence that non-newtonian flow properties have on the sliding. 


PACS. 87.19.rh Fluid transport and rheology - 83.60.Rs Shear rate-dependent structure - 47.50.Cd Non-newtonian 
fluid flows Modeling - 47.11.St Multi-scale methods 


1 Introduction mentally Il25ll26ll4ni42l and numerically/theoretically II43II44I 


Controlling and predicting the properties of spreading and mov¬ 
ing droplets is a major scientific challenge, relevant for an am¬ 
ple variety of applications, particularly in droplet-based mi¬ 
crofluidics |[T]|2]|3l4]|5]|6l|2llHl- From the point of view of funda¬ 
mental research ll9l[Tnl[TTl[ra[T3l[T4l[T^[T6l[T7lfT^^ . the spi 
ing and motion of a droplet on a solid surface has received 
much attention because of the singularity that occurs at the con¬ 
tact line El, where the competition between capillary forces 
and viscous dissipation results in a singular problem 015II16I 
El. Close to the contact line, the flow is confined to a wedge¬ 
like region whose shape is largely independent of the macro¬ 
scopic geometry. An added richness is brought by the complex 
non-newtonian behavior of the bulk phases, e.g. concentrated 
polymer solutions, which often exhibit both shear-thinning and 
normal stress effects 021ll22ll23ll24l . Deformation rates imposed 
on the fluid may indeed induce non-newtonian behavior within 
the wedge-like region 025II26I . Shear thinning was even pro¬ 
posed as mitigating the moving contact line singularity: if the 
viscosity of the liquid decreases monotonically with increas¬ 
ing shear stress, then there will still be an infinite shear stress, 
but the stress singularity will be integrable, yielding a finite 
force at the substrate Il27ll28ll29l . Some experimental studies 
analyzed the effects of non-newtonian rheology on the problem 
of spreading 0301131113211^ . The spreading behavior was found 
not to deviate strongly from that of the newtonian fluids E). 
The experiments were also compared with the predictions of 
lubrication theory of power-law fluids 035II36I . Other studies 
on non-newtonian effects focused on coating flows 037113811^ 
l40l and moving contact lines 025 ll26ll4ni42ll4Tl . both experi- 


l45l . Garoff and co-workers 0251126114411421 conducted a series 
of experiments on dynamic wetting: by comparing the steady 
interface shape, they observed that the non-newtonian effect 
was confined to the close vicinity of the contact line. In this 
inner region, the viscous bending of the interface is reduced by 
^ear thinning and enhanced by viscoelasticity, a fact that was 
also confirmed by numerical simulations mu. 

In this paper we report on a joint experimental/numerical study 
aimed to investigate the sliding of droplets made by polymeric 
solutions exhibiting non-newtonian properties. The correspond¬ 
ing newtonian problem has been studied in an ample variety 
of situations lfT11fT2lfT5lfT6lfT7lfT8lfT9ll20ll46ll47ll48l . but the ef¬ 
fects of viscoelasticity have been rarely investigated 049115011511 
I52ll5?ll54l . especially from the experimental point of view. As 
far as we know, only droplets of polystyrene solutions have 
been explored experimentally on hydrophilic substrates, where 
sliding is often affected by the presence of pearling ll55l . Most 
non-newtonian fluids exhibit various non-newtonian properties: 
notably, concentrated polymer solutions are characterized by 
both shear-thinning and normal stress effects 030113111^1^ . 
With the aim of concentrating on the shear-thinning flow prop¬ 
erties, we will consider solutions of Xanthan, a stiff rodlike 
polysaccharide. Xanthan solutions are chosen because, when 
compared to other flexible polymer solutions 03^ . are reported 
to have a strong shear rate dependence of the viscosity, and 
smaller elastic effects (i.e. normal stresses) 030113111^1^ . Op¬ 
eratively, the impact of a shear-dependent viscosity on the slid¬ 
ing process is analyzed through the relation between the av¬ 
erage dimensionless velocity (i.e. the Capillary number) and 
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the dimensionless volume forces (i.e. the Bond number). Re¬ 
sults from experiments are complemented with numerical sim¬ 
ulations of sliding droplets with the lattice Boltzmann models 
(LBM) 019II2OI . Numerical simulations provide an extremely 
versatile tool to investigate the effects of non-newtonian rhe¬ 
ology and to elucidate which deviations from the newtonian 
result could be attributed to a specific non-newtonian property. 
The paper is organized as follows: sectiondeals with the ex¬ 
perimental details, while the corresponding numerical results 
are described in section [3 Conclusions will follow in section 

SI 


2 Experiments 



We have investigated shear thinning solutions of Xanthan (molec¬ 
ular weight, M^^, ^ 10^ g mol^^ , Sigma Aldrich), a stiff rod¬ 
like polysaccharide. Solutions have concentrations ranging in 
the dilute or semi-dilute regime 032ll33ll56l as listed in Table 
[T] These non-newtonian solutions are mainly characterized by 
a shear-dependent viscosity, but also normal stress effects have 
been reported, particularly from high enough (~ 1000 ppm) 
concentrations . While stud 

ies on the shear viscosity of Xanthan gum solutions are well 
established, it seems from the literature that a smaller num¬ 
ber of studies have looked systematically and quantitatively 
at the first normal stress difference Ni ElSilEaiSllollSD 
I62ll63]l64ll65]l66ll^ . Specifically, some studies focused on the 
scaling of Ni vi. j for various polymer concentrations. Dif¬ 
ferent scaling laws were reported 057ll62IIF3ll64ll^ . depend¬ 
ing also on the range of shear rates used 057ll65ll66ll . In other 
works 1321, where the behavior of Xanthan solutions was com¬ 
pared to other polymeric solutions (e.g. polyacrylamide), nor¬ 
mal stresses were claimed to be negligible if compared to the 
strong shear rate dependence of the viscosity. 


CO 

Q. 


Liquid 

concentration 
(ppm w/w) 

a 

b 

Xanthan/water 

400 

0.033 

0.41 

Xanthan/water 

800 

0.076 

0.51 

Xanthan/water 

1500 

0.31 

0.61 

Xanthan/water 

2500 

0.99 

0.70 


described by the power law r]x, 
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^ Xanthan/water 500 ppm (Bonn, Meunier [59]) 

+ Xanthan/water 1000 ppm (Whitcomb, Macosko [57]) 
^ Xanthan/water 1500 ppm 

X Xanthan/(water+sucrose) 2000 ppm (Won, Kim [68]) 
O Xanthan/water 2500 ppm 
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Table 1: Viscosity parameters for the non-newtonian Xanthan 
solutions: a and b are the parameters fitting the rheological data 


n{j)=aj ° and plotted in Figure 


We characterized the rheological behavior of our solutions 
through a parallel plate rheometer (Advanced Rheometric Ex¬ 
pansion System ARES, TA instruments) equipped with two 50 
mm diameter disks. As shown in the top panel of Eigure[T] the 
viscosity varies with the shear rate 7 according to the power 
law 

rixanii) = 

with a and b dependent on the concentration, as listed in Ta¬ 
ble [T] (see also The bottom panel of Eigure[T] reports the 
evaluation of the first normal stress difference N\ achievable 


Eig. 1: Rheological measurements performed with a paral¬ 
lel plate rheometer. Top Panel: viscosity characterization for 
Xanthan solutions with different concentrations. The viscos¬ 
ity varies with the shear rate 7 according to the power law 
Tlxaniy) = 07 with a and b dependent on the concentra¬ 
tion, in agreement with 1^ . Bottom Panel: Ni measurement 
for 1500 and 2500 ppm Xanthan solutions in comparison with 
literature data 059II57II68II . 


by our rheometer, i.e. the values about the 1500 and 2500 ppm 
solution: values related to lower concentrations were found to 
be smaller, actually negligible below 1000 ppm l59l . and also 
a bit more scattered than data at higher concentrations because 
near to the sensitivity limit of the instrument. Our measure¬ 
ment indicates the presence of a measurable Ni in agreement 
with M57II68I . which investigate the rheology of Xanthan solu¬ 
tions of similar concentrations (see Pigure[T]i. 

The substrate used in all experiments was a homogeneous, poly¬ 
carbonate (PC) plate, whose wettability properties have been 
determined through the sessile droplet method. The sample is 
characterized by a static contact angle 9 close to 90° and dy¬ 
namic contact angles similar for all solutions, as reported in 
Table [2] 

Sliding measurements of 30 jJ-L droplets have been per¬ 
formed by depositing the desired volume on the already in¬ 
clined sample (see figure| 2 ]i and recording droplet motion with 
a CMOS camera (MV-D1024E Camera Link), whose maxi- 
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Liquid 

concentration 
(ppm w/w) 

0 

(degrees) 

Oa 

(degrees) 

(degrees) 

Ae 

(degrees) 

water 


81±3 

91±2 

66±3 

25±4 

Xanthan/water 

400 

86±2 

89±2 

64±2 

24±3 

Xanthan/water 

800 

87±2 

92±2 

67±3 

26±3 

Xanthan/water 

1500 

89±1 

91±1 

66±1 

25±2 

Xanthan/water 

2500 

87±3 

89±1 

63±2 

25±2 


Table 2: Static (9), advancing ( 0 a), receding (0^) contact angles and contact angle hysteresis (40) of all the investigated solutions 
on a homogeneous PC sample. 


mum frame rate is almost 120 fps 018111911^ . Acquired im¬ 
ages were analyzed through a custom-made program which 
first identifies the droplet contour and then fits it with a poly¬ 
nomial function, which is subsequently used to evaluate the 
front and rear contact points 069II7OI . After a transient time, 
the droplets acquire a stationary velocity that was evinced by 
the temporal evolution of the front contact point. For every so¬ 
lution, droplet sliding was observed at inclinations of the PC 
surface comprised between the sliding angle Uc (the smallest 
angle beyond which the droplet starts moving, see later) and 
the highest angle for which volume control was still accurate. 



Fig. 2: Sketch of the experiment: a droplet containing Xanthan 
at a given concentration is deposited on the polycarbonate sur¬ 
face, inclined by the angle a, through a vertical syringe pump 
dosing ~ 30 fiL with an accuracy of about 3%. The droplet is 
illuminated by a square 5x5 cm^ white LED backlight. The 
drop velocity U is measured by acquiring the whole sliding se¬ 
quence, g being the gravity acceleration. 


The main experimental results, consisting in the dependence 
of the stationary droplet velocity U on the inclination angle a 
of the PC surface, are shown in the top panel of figure |3] By 
comparing the results of the polymeric solutions with those of 
water, we observe a decrease of t/ as the polymer concentration 
increases due to the higher viscosity of the more concentrated 
solutions, which results in a smaller steady velocity at a fixed 
driving force 0111115 IflblfTTl . As a consequence, in order to bet¬ 


ter compare the behavior of the different liquids, the data have 
to be plotted in terms of dimensionless numbers which can be 
suggested by the analysis of the sliding of newtonian droplets 
mm, where the steady velocity U comes out from the bal¬ 
ance between the work done by the external driving force and 
the viscous dissipation. More quantitatively, a newtonian liquid 
droplet of volume V, density p and viscosity 77 sliding down 
an inclined plane tilted by an angle a is subject to the gravity 
force, the viscous drag and interfacial forces (see figure|2l). The 
down-plane component of the droplet weight is pgV sin a, the 
viscous drag force is proportional to VjV^I^U while the interfa¬ 
cial force is proportional to where a is the liquid-gas 

surface tension and Aq is a non dimensional factor depending 
on the perimeter-averaged projection of surface tension. In ad¬ 
dition, the non vanishing difference between the advancing and 
the receding contact angles for small velocities leads to contact 
angle hysteresis IfTTI . The critical inclination angle Uc is then 
identified as the angle below which the droplet stays pinned 
and it does not move. For inclination angles above Uc, the bal¬ 
ance of the gravity, viscous and capillary forces implies a linear 
scaling law II15II16I between the Capillary number Ca = TjU/(J 
and the Bond number Bo = g sin a/O’ 

Cao^Bo —Boc (1) 

where Boc = V^/^pg sin Ot/o depends on the wetting hystere¬ 
sis through Aq. 

For the non-newtonian Xanthan solutions the calculation of 
the Bond number is straightforward, on the other hand, due to 
the shear thinning behavior, the viscosity inside the droplet is 
obviously non constant and the computation of the Capillary 
number requires the evaluation of the viscosity corresponding 
to the sliding velocity U: one needs to devise a method to con¬ 
vert the data of viscosity versus shear-rate in viscosity versus 
U. For this practical purpose, we define an “effective” viscosity 
in terms of an effective shear rate Jeff = U/X 

= ajgff = a (2) 

where a and b are the parameters reported in Table [1] The def¬ 
inition of this effective viscosity hinges on the introduction of 
the characteristic scale X. The viscosity rjeff in Q is then used 
to calculate the Ca for the Xanthan droplets and the corre¬ 
sponding values are plotted as a function of Bo — Boc in the 
bottom graph of Figure [3] Remarkably, a constant value of X 
(~ 10 pm) for all the Xanthan concentrations investigated, al¬ 
lows to obtain a good collapse of the non-newtonian data on 






































4 


Silvia Varagnolo et al.: Sliding droplets of Xanthan solutions: a joint experimental and numerical study 




Bo - Bo 

c 

Fig. 3: Top panel: velocity of 30 /fL droplets sliding down the 
homogeneous plane (PC) inclined by different angles a for wa¬ 
ter and Xanthan/water solutions. Bottom panel: corresponding 
dimensionless Ca vi. Bo — Boc trend. 


the newtonian ones for small Bo. The scale X can be thought 
of as defining an effective gradient U/X which is “represen¬ 
tative” for viscous dissipation inside the droplet. Such scale 
is indeed smaller than the macroscopic length scale (i.e. the 
capillary length) of the order of 1 mm and larger than the mi¬ 
croscopic length (i.e. the slip length) ranging from 1 to 10^ nm 
ll43l . However, at larger Bo — Boc, non-newtonian data system¬ 
atically deviate from the reference newtonian data, and even a 
plateau is emerging for the largest concentrations considered. 
We think there may be different reasons for these deviations. 
First, the characteristic scale X may not be constant as a func¬ 
tion of the velocity; the non constancy of the viscosity inside 
the droplet may well induce a variation of X for larger veloc¬ 
ities, with X being an increasing function of the velocity. To 
shed some lights on this issue, we have carried out some nu¬ 
merical simulations with a power-law fluid and we have not 
found evidences of this velocity dependence in the character¬ 
istic scale A. These results are described in section 13.11 We 
could, instead, relate the deviations from non-newtonian data 
to the presence of normal stresses inside the droplet, which 


are more important at high shear rates 032ll57ll58ll59ll^l6ni62l 
I6^l^l65ll66ll67ll68l . i.e. at high velocities and Bo, while it is 
weaker at small inclinations where the behavior is similar to 
a newtonian fluid. With the measured velocities and the used 
values for A (10 /rm), we indeed obtain jeff ~ 10^ ~ 10^ s^* 
for the largest velocities of the more concentrated Xanthan so¬ 
lutions analyzed. This itself points to a representative value for 
the shear rate at which measured normal stresses are different 
from zero ll^l57ll5^l59ll6()ll^l62ll^lPll65ll^l67ll^ . as dis¬ 
cussed for the bottom panel of figure [T] An elastic non-viscous 
effect has been observed also in the flow of Xanthan solutions 
through a porous medium and has been ascribed to the pres¬ 
ence of normal stresses ll58l . 

To shed some lights on this issue, we have carried out some nu¬ 
merical simulations with a normal-stress fluid where thinning 
effects are extremely reduced. These results are described in 
section lT^ 


3 Numerical Simulations 

In our numerical study, we considered the diffuse-interface Navier- 
Stokes (NS) hydrodynamic equations of a binary mixture of 
two components 043ll71ll72ll73ll741175l obtained with a lattice 
Boltzmann model (LBM) lfT9ll to simulate a droplet (d) with 
dynamic viscosity rj^ and density sliding in an outer (o) fluid 
with dynamic viscosity rjg and density Po 1119112011 . We analyzed 
the case of a two-dimensional (cylindrical) droplet m with 
radius on a homogeneous surface with neutral wettability. 
The use of two-dimensional numerical simulations allowed us 
to better resolve the hydrodynamics inside the droplet. All the 
technical details for the simulation of newtonian droplets slid¬ 
ing down homogeneous surfaces are well described in our pre¬ 
vious publication im. A few words are in order for the choice 
of the viscous ratio % between the droplet and outer phases: ap¬ 
proaching the extremely large experimental liquid-air viscous 
ratio is not easily achievable with the numerical models. One 
can partially cope with this problem by making the body force 
monotonously decrease to zero across the fluid-fluid interface, 
accounting for the fact that the outer phase remains inert in the 
limit of zero droplet size ll20l . In fact, when working with very 
large viscosity ratios, numerical errors increase close to the in¬ 
terface ||43]| . especially close to the contact line. Fortunately, 
these deficiencies do not affect conclusions with viscous ra¬ 
tios II19II43I up to a few tens. 

3.1 Shear Thinning Fluid 

The numerical method, based on LBM, offers a unique advan¬ 
tage to attack these problems, being able to integrate simulta¬ 
neously the non-trivial interface dynamics (due to the diffuse 
interface nature of the method) and the non-newtonian nature 
of the bulk phases, by making the relaxation times in LBM de¬ 
pendent on the local shear properties Ct). The equations we 
solve in the outer (o) and droplet (d) phases are the NS equa¬ 
tions 

Po,d T {V‘o,d ' = ^Po,dX- 
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Here, d and rjo ii are the velocity and the dynamic viscosity 
of the outer and droplet phase, respectively. The correspond¬ 
ing densities and pressures are indicated with Poj and 
while {'^UociY indicates the transpose of (Vmo,(/). As plot¬ 
ted in Figure 01 the viscosity in the droplet phase is chosen to 
be 


Idir) 


i?o + fl7o 7<7o 

rio + af-‘’ 7>7o 


(4) 


where a = 0.01 (Ibu, lattice Boltzmann units) and where the 
cut-off shear rate 70 = 10 ^^ Ibu is introduced to prevent a di¬ 
vergence of the viscosity at small shear rates. The background 
viscosity rjo = 0.07 is chosen to avoid a zero viscosity in the 
high shear rate limit, which would cause numerical instabili¬ 
ties. For all practical purposes, in the range of shear rates rele¬ 
vant for our problem (see figure 01 i, we are actually simulating 
a power-law fluid. In the outer matrix phase the viscosity is 
kept constant to rjo = 0.07 so that the viscous ratio between the 
droplet and the outer phase ranges in the interval [1 — 10 ] for 
the results presented. 






10 


- b=0.0 newtonian 

-b=0.2 

. b=0.4 

- b=0.6 

1 

^ 1 

0 

0 

w 

■> 

ro 0.1 


0) 

w 

0.01 





10 " 10 “ 10 " 10 " 10 “ 
shear rate (Ibu) 


Fig. 4: Viscosity of the power-law fluid used in the numerical 
simulations. The exponent b regulates the degree of thinning. 
On average, dissipation is larger at increasing b, in qualitative 
agreement with the behavior of the viscosity of Xanthan solu¬ 
tions at changing the polymer concentration. 


In the top panel of figure |5] we report the stationary veloc¬ 
ity measurement for 3 cases: a newtonian case {b — 0) and two 
viscoelastic cases with b — 0.2 and b = 0.4. Consistently with 
the fact that the viscosity chosen (01 produces a more viscous 
fluid at increasing b, the stationary droplet velocity is smaller 
at increasing the thinning exponent. However, by following the 
same procedure explained for the experimental data, once we 
tune the scale X to achieve good collapse on the newtonian 
data for small velocities, a very good collapse is obtained also 
for higher Ca, as shown in the bottom graph of Figure |5] This 
surely indicates robustness in this “practical” procedure. It is 
true that some tiny deviations emerge for the largest Ca, but 
they are really very small and surely not at all reminiscent of the 
plateau behavior observed for the experimental data on Xan¬ 
than. 
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Fig. 5; Top panel; stationary velocity for sliding droplets at 
changing the thinning exponent b in (01i. Bottom panel: rescaled 
data based on the lengthscale X = 0.3 (Ibu) in the definition of 
the shear rate j=U/X used in the viscosity to identify the “ef¬ 
fective” viscosity of the thinning fluid. The rescaling is chosen 
in such a way that the non-newtonian data overlap the newto¬ 
nian data for small Bond numbers. Good collapse is also found 
at larger Bond number. 


3.2 Normal Stress Fluid 

In order to introduce normal stresses in the droplet phase, we 
used an available version of our code implementing the FENE- 
P constitutive model, obtained via a pre-averaging approxima¬ 
tion applied to a suspension of non interacting Einitely Extensi¬ 
ble Nonlinear Elastic (EENE) dumbbells. This is well-adapted 
for dilute (and semidilute) polymeric solutions and was pre¬ 
viously used to analyze filament thinning of polymeric fluids 
in macroscopic experiments II23II78I and droplet breakup pro¬ 
cesses in microchannels El. EENE-P fluids exhibit normal 
stress effects with a quadratic dependence on the shear rate. 
Such dependence on the shear rate contrasts with the sub-linear 
scaling law found for rigid polysaccharide solutions like Xan¬ 
than 162 ; in fact, we rather use this numerical model to “gener- 
ically” highlight the effect of Normal Stresses. The algorithm 
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we use is based on a hybrid combination of LBM with finite 
difference schemes, as detailed in II79II80I . The equations that 
we solve in the droplet phase are the Navier-Stokes (NS) equa¬ 
tions coupled to the FENE-P constitutive equations 

Pd [dfUd + (ud ■ V)Mrf] = -VPd+ 

V {Pd{^ud + (yudf)) + ^V • [f{rp)C]. 

+{ud ■ V)C = c ■ {Vud)+{^Udf ■ c 

_ l^ f{rp)C-I \^ (6) 

where rjp is the viscosity parameter for the EENE-P solute, 
Zp the polymer relaxation time, C the polymer-conformation 
tensor, I the identity tensor, f(rp) = (L^ — 3)/(L^ — rj,) the 
EENE-P potential that ensures finite extensibility, rp = y^Tr(C) 
and L is the maximum possible extension of the polymers ED 
l22l . In the outer phase we just consider the NS equations with 
viscosity r\p -b 77 ^, i.e. with the same shear viscosity as the 
droplet phase. The chosen value of is so large that the sys¬ 
tem only displays very weak thinning effects, while possessing 
normal stresses increasing with the squared shear rate (see fig- 
ure|6]l. The value of the polymer relaxation time Zp is chosen in 
such a way that the range of sliding velocity produces shear rate 
in the wedge which are of the order of magnitude of Zp^ . We 
expect that this choice allows us to appreciate the effects of the 
normal stresses on the Capillary-Bond curve. Eurther details on 
the algorithm used can be found in 17911801 . 



10 '' 10 " 10 ^ 10 ’ 10 “ 10 ’ 
shear rate (Ibu) 


Eig. 6 : Normal Stresses (main panel) and shear viscosity (inset) 
curves for the chosen viscoelastic solution. With the parame¬ 
ters chosen, the thinning effects are minimized while normal 
stresses are present for the range of shear rates relevant for our 
problem. 


Results for the Ca vi. Bo curve are presented in figure |7l 
The concentration Pp/Ptot — Pp!{P d + Ip) N systematically 
varied from 0.01 to 0.15. The viscosity ratio of the two fluids 
was chosen in such a way that x = (Vd + Tjp) jrjo = ^ -0. This 
allowed to study the properties of viscoelastic droplets (rjp ^ 0 ) 


and compare them with the corresponding newtonian (rip = 0 ) 
case. It is readily verified that the deviations from linearity first 
emerge with a slight sublinear behavior for the smallest con¬ 
centrations, while they produce a plateau for the largest con¬ 
centrations. The behavior observed in the presence of normal 
stresses points to the fact that effects other than viscous dis¬ 
sipation come into play to balance the work done by the ex¬ 
ternal driving. We have verified in the numerical simulations 
that normal stresses provide an additional driving and accel¬ 
erate the flow near the contact line ll^ . and the effect does 
so more strongly as the shear rate is increased. Parts closest 
to the rear wedge move faster and there the polymer feedback 
stress is enhanced to provide a counter force against sliding. It 
is interesting to observe that the characteristic Ca at which the 
linear behavior starts to be violated is essentially unchanged. 
Preliminary numerical simulations show that such “critical” Ca 
depends on the relaxation time Zp: numerical results obtained 
at changing the relaxation time Zp will be the object of a forth¬ 
coming publication, which will be useful to better characterize 
the relative importance of both the polymer relaxation time and 
the polymer concentration in generating sublinear and plateau 
behaviours. 


0.02 


0.01 


0.00 


Pig. 7: Capillary vw. Bond number curve for viscoelastic fluid 
with rheological properties described in figure | 6 l The concen¬ 
tration rjp/rjpoT = ’np/i'Hd + Pp) N varied from 0.01 to 0.14. 
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4 Conclusions 

We have investigated the sliding of droplets made of polysac¬ 
charides in aqueous solutions. In particular, we used Xanthan, 
a stiff rodlike polysaccharide exhibiting a non-newtonian be¬ 
havior, notably characterized by a shear-rate dependence of the 
viscosity. The sliding of non-newtonian droplets was compared 
to the one of water through the relation between the Capillary 
number (i.e. the dimensionless velocity) and the Bo number 
(i.e. the dimensionless volume force), as the driving force is 
changed by varying the plane inclination. To account for the 
spatially varying viscosity inside the droplet, we argued that an 
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effective viscosity could be defined to quantitatively compare 
with the corresponding newtonian problem. The usual linear 
behavior of newtonian solutions is found to well adapt to the 
non-newtonian solutions only up to a given Bo number, above 
which deviations emerge, as the linear behavior is first replaced 
by a sublinear behavior (small concentrations) and then by a 
more visible plateau exhibiting a saturation in the Capillary 
number. To discriminate to what degree the observed behav¬ 
ior could be due to the thinning effects or other non-newtonian 
effects (e.g. Normal Stresses), results from experiments were 
complemented with lattice Boltzmann numerical simulations 
of sliding droplets, aimed to disentangle the influence that non- 
newtonian flow properties have on the sliding. Numerical sim¬ 
ulations reveal that the observed behavior is more likely at¬ 
tributed to the emergence of normal stresses inside the non- 
newtonian droplet, rather than thinning effects. 

For future studies it is definitively warranted a systematic in¬ 
vestigation of the sliding problem for other polymeric solu¬ 
tions. This is the case, for example, of polyacrylamide solutions 
where the viscosity is known to be only weakly dependent 
on the shear rate while being characterized by normal stress ef¬ 
fects. Work in this direction is ongoing. 

The authors kindly acknowledge funding from the Euro¬ 
pean Research Council under the European Community’s Sev¬ 
enth Eramework Programme (EP7/2007-2013)/ERC Grant Agree¬ 
ment No. 279004. We are particularly grateful to Prof. Gio¬ 
vanni Lucchetta (Department of Industrial Engineering, Uni¬ 
versity of Padova, Italy) for fruitful discussions and support 
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